Scaling Distributions of Quarks, Mesons and Proton for all px, Energy and Centrality 
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We present the evidences for the existence of a universal scaling behavior of the production of tv° 
at all transverse momenta in heavy-ion collisions at all centralities and all collision energies. The 
corresponding scaling behavior of the quarks is then derived just before the quarks recombine with 
antiquarks to form the pions. The degradation effect of the dense medium on the quark pr is derived 
from the scaling distribution. In the recombination model it is then possible to calculate the pr 
distributions of the produced proton and kaon, which are scaling also. Experimentally verifiable 
predictions are made. Implications of the existence of the scaling behavior are discussed. 

PACS numbers: 25.75.Dw, 24.85. +p 



I. INTRODUCTION 

In two recent papers we have discussed the scaling 
properties of the large pt distributions of 7T° produced 
in Au+Au collisions at the relativistic heavy-ion collider 
(RHIC) and presented their implications. As reported in 
the first paper Q], hereafter referred to as I, we found 
energy scaling at maximum centrality, while in the sec- 
ond paper 0, referred to as II, we found centrality scal- 
ing at the highest energy. The two can be combined to 
yield one scaling distribution for all energy and centrality. 
In I we derive the quark distribution from the 7r° data 
in the framework of the recombination model; we now 
do the same for all centralities and determine the na- 
ture of degradation of the quark momentum in the dense 
medium, as we have done for the 7r° momentum in II. 
From the quark distribution we can calculate the proton 
distribution for all centralities. Moreover, we can extend 
our consideration to the production of kaons so that we 
can calculate not only the p/ir ratio, but also the K/w 
ratio. 

We are able to do all that for two essential reasons. 
The first is that the discovery of the scaling properties 
facilitates the analysis by avoiding the need to consider 
the variation of physics issues at different centralities and 
energies. The other is that the quark distribution we de- 
rive is for q and q just before hadronization. It is not the 
result of some dynamical evolution starting from hard 
collisions, for which many complex issues must be con- 
sidered y, El EJ . The recombination model that we use 
can only address the hadronization problem of the soft 
partons at low virtuality, but at any px- From the pion 
data we infer the distributions of the soft q and q, which 
in turn are used to give the proton and kaon distribu- 
tions through recombination. How the soft partons get 
to be where they are in the pt space is not considered. 
However, by studying the momentum degradation of the 
quarks, we gain from the centrality dependence of the px 
distributions some understanding about how quarks lose 
momenta as they propagate through the dense medium. 
Experimentally verifiable predictions are made on the 



proton and kaon transverse momentum distributions. 

The physical interpretation of the scaling variable is 
given at the very end, where the term transversity is 
suggested to refer to the difficulty of acquiring trans- 
verse motion. The broader implication on the creation of 
quark-gluon plasma is finally addressed. 



II. SCALING DISTRIBUTION OF PIONS 

From the preliminary PHENIX data of ir° production 
in Au+Au collisions at RHIC we have found a scaling 
distribution at midrapidity 



$0) = A(N)K 2 (s, N) 



1 cPNk 
2ttpt drjdpT 



where 



z = p T /K{s,N) 



(1) 



(2) 



The symbol N denotes the number of participants, A part , 
for brevity. The scaling factor K(s, N) is first found in 
I for N fixed at its maximum, JV max = 350, to be (in 
GeV/c) 



K(s) = 0.69+ 1.55 x 10~ 3 \/i , 



(3) 



where yfs is in units of GeV, and K(s) is normalized to 
1 at y/s = 200 GeV. When N is allowed to vary, while 
y/s is fixed at 200 GeV, it is found in II that the scaling 
factor is 



K(N) = 1.226 - 6.36 x 10~ 4 A 



(4) 



normalized to 1 at N — 350. We now combine the two 
and assume the factorizable form 



K(s,N) = K(s)K(N) 



(5) 



In this paper we investigate the centrality dependence 
mostly at y/s = 200 GeV. The normalization factor A(N) 



2 



in Eq. QJ is found in II to have a power-law dependence 
on the number of binary collisions, N c , 

A{N C ) = 530N- 9 , (6) 

where N c in turn depends on N as 

N c = 0.447V 1 - 33 , (7) 

a relationship that is determined from the tables listed 
in Refs. 
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FIG. 1: Scaling distribution of pion production. Data at 
y/s = 130 and 200 GeV are preliminary and were presented 
by PHENIX 0, for Au+Au collisions. The data at 17 
GeV are from WA98 for Pb+Pb collisions. The solid line 
is a fit by Eq. ©. 

In Fig. 1 we show the combined plot of $(z) exhibiting 
the y/s = 200 GeV PHENIX data for 5 bins of centrality 
as well as the y/s = 130 GeV data at 0-10% centrality [f| 
and the 17 GeV data at 0-12.7% 9]. The tt° data, which 
are only for pt > 1 GeV/c, are supplemented by the 7r + 
data forpT < 1 GeV/c at y/1 = 200 GeV and 0-5% cen- 
trality 10]. Evidently, all the data points fall on one uni- 
versal scaling distribution that is invariant under changes 
in N and y/s. The 17 GeV data are obtained by the 
WA98 Collaboration 9] for Pb+Pb collisions and show 
a slight departure from the universal curve for z > 3. It 
should be recognized that those data points that devi- 
ate from the scaling distribution correspond to pt > 3 
GeV/c, which is a pr range that represents a very large 
fraction of the total available energy of 17 GeV. Thus the 
kinematic constraint of energy conservation introduces a 
non-dynamical factor that suppresses the high-p-r behav- 
ior, not present in the other data at y/s >130 GeV. Such 
a violation of scaling is expected, and should not be re- 
garded as an invalidation of the general scaling behavior 
we observe. On the contrary, it is amazing that the scal- 
ing behavior can cover such a wide range of y/s, when 
most of the 17 GeV data points with pr < 3 GeV/c are 
included. 

The scaling data points can be well fitted by 



which is shown by the solid line in Fig. 1. This formula 
differs from the one given in I mainly by the addition 
of the exponential term, which reflects the statistical be- 
havior at small z represented by the tt + data, but is in- 
significant at large z, where the power-law behavior is 
indicative of the effects of hard collisions. 

A number of consequences of the scaling distribution, 
Eq. JHJ, can be examined directly. First, the integral 



2 = 



dzz^(z) 



A(N) dN^o 
2n drj 



(9) 



can be evaluated to yield I = 46.2. Using Eqs. ijfjj and 
for A(N) = A(N C (N)), we can calculate the N de- 
pendence of dN^o /drj at midrapidity. The data that are 
available for comparison are dN c h/dri/(0.5N) vs N at 
y/s = 130 GeV 0, shown in Fig. 2. We plot in that figure 
in solid line our calculated result for dN^± jdr\j (0.5N), 
which is obtained from Eq. © by identifying dN n ± /drj 
with 2dN n o/dr]. Since tt^ is the dominant part of all 
charged particles, one should regard the comparison be- 
tween the calculated N^± and the measured iV c h to be 
satisfactory. Recall that although the p/ir ratio can ex- 
ceed 1 around px ~ 3 GeV/c, it is small at small pr 
where the distributions are dominant, so the production 
of proton does not contribute to the integrated result as 
a large fraction of the total dN^/drj. 
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FIG. 2: Calculated dN n ± /dr)/(0.5N) compared to the data 
on dN ch /dri/(0.5N) @. 

Another consequence of &(z) is the possibility to cal- 
culate the nuclear modification factor 



Raa(pt) 



d'NX/dydpr 
N c d 2 N™/dr)dp T 



(10) 



$(z) = 1200 (z 2 + 2) 



-4.8 



l + 25e- 4 ' 5z ) 



(8) 



The data available for that are given in Ref. p. We 
show in Fig. 3 the data for two centrality bins: 0-10% 
and 70-80%. The corresponding values of N are 317 and 
9.5, respectively 0, for which the calculated Raa(pt) 
are shown by the solid and dash-dot lines. For N p f, in 



3 



the denominator of Eq. 1|10[) we have used N — 2. Ex- 
perimentally, it is known that the peripheral nuclear col- 
lisions cannot be identified with pp collisions. The ratio 
Raa{pt) is defined with N p „ in the denominator so as 
to have a definitive experimental normalization. We cal- 
culate the ratio with N = 2 in the denominator so that 
no additional experimental input is used. Our point is 
to show the consistency of $(z) when the value of N is 
extrapolated to extreme limits. Toward that end we find 
the comparison between the data points and the calcu- 
lated curves in Fig. 3 to be acceptable. 



< 
< 
DC 



1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 

0. 



N=9.5 
N=317 



•70-80% 
°0-10% 








4 



FIG. 3: Nuclear modification factor Raa(pt) for 7r° produc- 
tion. The preliminary data are from PHENIX ;G|. Fhe lines 
are calculated results using N = 2 extrapolation of 4>(z) for 
normalization rather than using independent pp distribution. 

The scaling distribution $(z) is a representation of the 
data over 9 orders of magnitude for all centralities and 
for all energies where data exist. Since the fit is done in 
the log scale in Fig. 1, one can expect some deviations 
in the linear scale, as the tests in Figs. 2 and 3 are 
done. The data must, of course, be self-consistent, so any 
discrepencies in those figures are due to the extrapolation 
of $>(z) to very small z in the calculation of dN^a/dtq and 
to very small N for Raa- We conclude from those tests 
that the formula (JSJ) for Q(z) is quite reliable even down 
to very small values of z and N. 



III. SCALING DISTRIBUTION OF QUARKS 

From the scaling distribution of 7r° it is possible to 
derive the quark distribution in the framework of the re- 
combination model. Those quarks have low virtuality 
(hence, soft) and are at the last stage of their existence 
just before hadronization. They are not to be confused 
with the partons at high virtuality (hence, hard) just af- 
ter hard collisions. The evolution from the hard partons 
to the soft quarks through gluon radiation and conver- 
sion to quark pairs in the dense medium involve both 
perturbative and non-perturbative QCD processes that 



are complicated, only some of which can be calculated 
H H GJ, E3 The last step to hadrons is circumvented 
by use of the phcnomcnological fragmentation function 
that connects hard partons to hadrons directly. Our use 
of the recombination model treats only the last step from 
soft quarks to hadrons, and can make no statement about 
the evolutionary process that begins from hard collisions. 
What can be treated is how the soft quarks recombine in 
different combinations to form different hadrons. That is 
what we shall do in the following sections. Here we first 
derive the quark distributions and examine how they de- 
pend on centrality. 

The application of the recombination model 0, 0] to 
the high px problem has been discussed in I. The recom- 
bination of a q and a q to form a 7r°, where qq can be 
either uu or dd, but not both uu and dd, is described by 
[see I-Eq. 0] 



diNLo 



" f 

j— = I dz x dz 2 z 1 z 2 F qq -(z 1 , z 2 )R„{zi, z 2 , z) , (11) 



zdz 



where our dN„o / zdz is averaged over rapidity at midra- 
pidity and over the azimuthal angle <p with the 1/27T 
factor included, unlike the experimental distribution 
{2ir)~ 1 d 2 N^o /dripTdpT, which is an average over tfi that 
shows the 1/2-k factor explicitly, as in Eq. QJ. The 
joint distribution F q q(z±, z 2 ) is assumed, for heavy-ion 
collisions, to have the factorizable form F q q[z\,z 2 ) = 
F q (zi)Fq(z 2 ) with Fq(z\) being the quark distribution 
in the scaling variable Z\, and Fq(z 2 ) for the antiquark. 
The recombination function R^{zi, z 2 , z) depends on the 
wave function of the constituent quarks in the pion; in 
the valon model Q, 0] it is 



Rtt{z\,Z 2 ,z) 



fz_ 1+ z 2 

\ z z 



- 1 



(12) 



where the valon distribution in the pion is determined by 
use of the data on Drell-Yan production by pion, which 
is the only way to probe the pion structure. 

At N = iVmax, considered in I, the LHS of Eq. JTTJ 
is identified with $(z), which in I is scaling in s. Now 
for all centrality the LHS of Eq. (|11|) is replaced by the 
new <f>(z) given in Eq. (JSJ), scaling in both N and s, and 
on the RHS the new scaling F q (zi) and Fq(z 2 ) are to be 
determined. Thus putting the various pieces together, 
we have 



d>(z) 



dz^fl-^Fqiz^Fqiz-zJ . (13) 



There is no explicit dependence on N or s in this equa- 
tion, but the q and q distributions, being scale invariant, 
have implicit dependences on pt and N that will be ex- 
amined below. 

Since the q and q in Eq. ifHfl) are at the end of their evo- 
lutionary processes, and are therefore soft partons domi- 
nated by the products of gluon conversion, their distribu- 
tions can differ in normalization, but not significantly in 
their z\ and z 2 dependences, a property that is supported 



4 



by the observation that the p/p ratio is nearly constant 
in pt fl5| . Denoting the p/p ratio by c, we thus use the 
relationship 



F q (z) = c^F q {z) 



(14) 



At RHIC energies c is roughly 0.7; at SPS c is about 0.13 
[lH Il7j . Since c depends on the energy y/s, Eq. lfTI)l is 
not strictly a scaling relationship. The smallness of c at 
SPS energy is due to the difficulty of producing a large 
number of p compared to p at y/s — 17 GeV 18]. Such 
a scaling violation is expected, just like the deviation of 
the diamond points at z > 3 from the universal scaling 
curve in Fig. 1 for kinematical reasons. However, over the 
whole range of variation of F q (z) that we shall determine, 
the effect due to the variation of c is small by comparison, 
as we shall see. 
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FIG. 4: Comparison of fitted and calculated curves of &(z). 
The solid line is a plot of Eq. |0J obtained from the fit in Fig. 
1. The dashed line is the result of the calculation using Eq. 

Let us put our emphasis on the scaling region of c by 
setting c = 0.7. Using that in Eq. (|14f> and then on the 
RHS of i fEfy . with Eq. JSJl on the LHS, we can vary the 
parametrization of F q (z) to achieve a good fit of §>(z). 
Our result is 



F q {z) = 15(z 2 + 0.47z + 0.72) 



-4.25 



(15) 



In Fig. 4 we show a plot of Eq. (JSJ by the solid line 
and our fit of it using Eq. (|15fl in (|13fl by the dashed line. 
Clearly, the fit is very good. The quark distribution given 
by Eq. I|15fl is shown in Fig. 5 over a range of 9 orders 
of magnitude of variation. How would that be affected, 
if c is lowered to 0.13? The normalization of F q (z) in 
Fig. 5 would be increased by a factor of (0.7/0. 13) 1 / 6 
I ..'!. which is insignificant compared to the 9 orders of 
magnitude of variation of F q (z). 

From Eq. (|15|l we can first determine the average z 



J dzz 2 F q (z) 
J dzzF q {z) 



= 0.415. 



(16) 
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FIG. 5: Scaling distribution of quarks before hadronization. 

We then define a new variable 

u = z/(z) , (17) 
in terms of which we can define a new distribution 

V>,(u) = F q {z{u)) / / duuF q {z{u)) (18) 



(z) 2 F q {z{u)) / / dzzF q {z) 



This distribution has not only the property that 



du uip q (u) 



by definition, but also 



du u 2 ipq{ u ) = 1 



(19) 



(20) 



These are the properties of a KNO-type distribution . 
To be strictly KNO scaling, all higher moments should be 
independent of N, as we shall investigate. The virtue of 
the scaling variable u is that it can be expressed directly 
in terms of px, since from Eqs. J2Jl and (|17fl we have 



: prf (p 



'TIN 



(21) 



where (pt) n is the average pt at fixed centrality and s. 
Note that the scale factor K(s, N) is common for the pt 
of both 7T° and q. In the case of pions [2| the variable u 
can, in principle, be determined unambiguously from the 
experiments directly, unlike the variable z that requires 
rescaling and fitting of the data at each N and s. Indeed, 
the variable u is constructed in the same spirit as in the 
original derivation of KNO scaling 0] . Although in the 
case of quarks here u cannot be experimentally measured, 
it is useful to have a KNO distribution ipqi u ) as a § oa l 
for theoretical modeling, since the u variable can more 
directly be related to the px of the quarks. 
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FIG. 6: KNO-type distribution of quarks. 
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FIG. 7: Centrality dependence of normalized quark distribu- 
tion at yfs = 200 GeV. 



From the definition in Eq. (18) we have 

tp q (u) = 1420 (u 2 + 1.13u + 4.18)" 4 - 25 . (22) 

A plot of this distribution is shown in Fig. 6. The exis- 
tence of such a scaling distribution in terms of an intu- 
itive variable given in Eq. I|21l) suggests that there is a 
great deal of regularity in the interplay between iV and 
Pt- Remembering how J/ip production is expected to 
have an anomalous px dependence in AA collisions when 
deconfinement occurs |20j , we find the lack of irregularity 
in the pt distribution of the quarks to suggest that the 
subject of pt dependences of low-mass particles is not a 
fertile ground to find signals of quark-gluon plasma — 
unless, of course, scaling violation is dramatically found 
at LHC, say. 

Having analytic forms for F q (z) and ij) q (u) enables us 
to investigate the degradation of the parton px in the 
dense medium. To that end we define the momentum 
fraction variable 



where Kq is a fixed scale, which we take to be 10 Gev/c. 
It is tacitly assumed that there is no physics of interest 
here for pt > 10 GeV/c. If that is not the case, an 
upward revision of Kq is trivial. From Eqs. J2J and l)23|) 
we have 

z = xK /K(s,N) , (24) 

and we may rewrite the scaling quark distribution as 

F q (z) = F q (x,N) , (25) 

where the s dependence is suppressed. We can then de- 
fine the normalized quark distribution 



Q(x,N) =F q (x,N) 



dx xF q {x, N) 



(26) 



The x dependences of Q(x, N) for various representative 
values of N are shown in Fig. 7 for y/s = 200 GeV. It is 
evident that, as N increases, the high x tail of Q(x,N) 
is suppressed, with a concurrent slight increase at very 
small x, since the integral of Q(x, N) over xdx is constant 
at 1. The crossover is at x ~ 0.04, corresponding to 
Pt ~ 0.4 Gev/c. Thus the effect of the dense medium is 
to degrade the p^ of the quarks, which is a well known 
property, but in a very regular way that we now describe. 
Let us define the ratio 



R(x,N) 



Q(x,N) 
Q(x,2) 



(27) 



so that the suppression at most x and enhancement at 
small x can be exhibited more clearly, as shown in Fig. 
8 for v/i = 200 GeV. Note that for x > 0.3 the suppres- 
sion is rather uniform. The rapid change in the range 
< x < 0.3 can be seen in a different plot, shown in 
Fig. 9. The decrease for x > 0.04, as TV is increased, is 
now clear, as is the increase for x < 0.04. It should be 
recognized that in normalizing Q(x, N) by Q(x, 2) in Eq. 
(|27|l it is not important whether Q(x, 2) agrees well with 
the corresponding distribution in pp collisions. The ra- 
tio removes the exponential dependence at small x that is 
common for all N and displays better the relative change 
as N is varied. Also, the extrapolation to very small x 
may not be accurate. The line for x — 0.01 in Fig. 9 
is intended mainly to give a rough idea of the nature of 
increase for x < 0.04. 

To quantify the degree of degradation, we take the mo- 
ments of Q(x,N) and define 



n (N) = / dxx n+i Q(x,N) . 



(28) 



x = pt/K 



(23) 



We do not consider the order of the moment n > 5, 
since higher moments demand more accuracy at higher x, 
which we cut off at x = 1 . The extrapolation of our scal- 
ing distribution F q (z) to higher z is not without uncer- 
tainties. However, for n < 5 our analysis is reliable, and 
provides adequate insight into the nature of the degra- 
dation in the < x < 1 region. In Fig. 10 we show 
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FIG. 8: Ratio of quark distribution relative to N = 2, ex- 
hibiting the degradation at x > 0.04 and the enhancement at 
x < 0.04 for ^fs = 200 GeV. 




Since Eq. (gHJ implies Q n (N) 
fore also 



IN 



, Eq. is there- 



-Xn(N-No) 



{x'-) N = {x") No e ""■ v " . (33) 

In particular, for n = 1, iVo = 2, and A*" = 350, we have 

{x) 350 /(x) 2 = 0.83 . (34) 

Thus even at the maximum separation between N and 
No, the average px of the quarks loses only 17%. This is 
because even severe suppression of high- a; quarks cannot 
change significantly the average (x) , which is dominated 
by the low- a; behavior. However, the same cannot be said 
about the higher moments. 

What we have found above are properties of the quarks 
before hadronization. Unfortunately, they cannot be 
checked directly by experiments. For testable predictions 
we now go to the study of proton and kaon formation, 
whose spectra can be measured experimentally. 
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FIG. 9: Same as Fig. 8 but plotted for various fixed values of 
x. 



lnQ n (N) vs N for n = 1, • • • , 5; the relationship can be 
well approximated by linear dependence 



lnQ„(JV) = a n -b„N , 



(29) 



where the slope b n is shown in the inset. Clearly, b n 
depends linearly on n 



b n = Xn, A = 5.35xl0^ 4 . 
We can combine these two equations to write 
d 



dN 



In Q n (N) = -\n 



n < 5, 



or 



Qn(N) = Q n (N )e- Xn ^ N - N ^ . 



(30) 



(31) 



(32) 



FIG. 10: Centrality dependence of the moments, Q n (N), 
whose log values are raised by the quantities in the paren- 
theses. The inset shows the slopes b n , the line being a linear 
fit. 



IV. SCALING DISTRIBUTION OF PROTON 

The quarks considered in the preceding section can not 
only combine with antiquarks to form pions, but also 
combine with other quarks to form protons. The formu- 
lation of the problem is discussed in I, where centrality 
variation is not considered. Now we allow N to vary, 
starting from the new quark distribution F q (z), whose 
scaling behavior includes the implicit dependence on cen- 
trality. As before, the treatment of proton production 
is not reliable at low z, where pt of the proton is low 
enough to make the mass effect important. Our calcula- 
tion that is scale invariant cannot take into account the 
mass-dependent effects. 
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The proton distribution arising from the recombination 
of uud quarks is given by [l| 



dN, 



zdz 
where 



and 



dz x dz 2 dz 3 z 1 z 2 z 3 F uud (z ll z 2 , z 3 )R p (zi, z 2l z 3 , z) , 

(35) 



F uu d(zi,z 2l z 3 ) = F u (zi)F u (z 2 )F d (z3) , (36) 



R p (z 1 ,z 2 ,Z3,z) = 0.057^-^(^1,6,6) ■ (37) 



Gp(6,6;6) is the valon distribution in a proton, ex- 
pressed in terms of the valon momentum fractions £j = 
Zi/z |2l| |. It is determined from the parton distributions 
that fit the deep inelastic scattering data, and is 

G P (6,6,6)=5(66) Q ^(£i+6 + 6-i) , (38) 

where 

a = 1.755, = 1.05 , (39) 



g = [B(a + l,/3 + l)B(a + 1, a + j3 + 2)}~ 



(40) 



The recombination function for the proton is more com- 
plicated than that for the pion because the proton is not 
as tightly bound as the pion in terms of the constituent 
quarks masses, but the procedures for the determination 
of the recombination functions are similar. 

Equation i|35l) is derived in I for N = N miix . For N < 
]V m ax the new scaling function for F q in Eq. (jT5|l is to be 
used for the quark distributions F u and Fd in Eq. (|3*H|l . 
The resulting integral is to be identified with the new 
proton scaling distribution <5> p (z), as in the pion case. 
We show in Fig. 11 the result on & p (z), for which only 
the z > 3 part is plotted, since it is unreliable for z < 3 
due to the neglect of the mass effect. 
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The relationship between the calculated & p (z) and the 
Pt distribution dN p /pxdpx{s,N) is 



%(z)=A p (N)K 2 p (s,N) 



dN p 
p T dp: 



-(s,N) 



(41) 



where, we repeat, our dN p /pxdpT corresponds to the ex- 
perimental (2-kpt)^ 1 d 2 N p / drjdpT [cf. Eq. QJ]. We have 
added a subscript p to each function in Eq. (|4*T)l to em- 
phasize their reference to proton; indeed, we should sim- 
ilarly add a subscript n to the corresponding functions in 
Eq. to clarify their differences, as we shall do below 
when we compare p and tt production. Whereas Eq. Q 
is determined by the phenomenological analysis of the 7r° 
data in II, there are no similar data on the p spectra to 
confirm Eq. (|41|) . However, on theoretical grounds we ex- 
pect that relationship to exist for the following reasons. 
First, since the quark distribution is scale invariant, the 
recombination model implies that there is a scaling dis- 
tribution &p(z) for the proton, as we have calculated. 
That gives the LHS of Eq. ijHJ). On the RHS we ex- 
pect K p (s,N) = Kt T (s,N) = K q (s,N), since the same 
scaling variable z has been used for pion, quark and pro- 
ton. Without that universal variable z the recombination 
model cannot be formulated in the form of Eqs. (|llfl and 
(|35|l . Finally, we conjecture that 



A p (N)=All\N) 



(42) 



on the grounds of internal consistency, since the com- 
parison among Eqs. (Q, (fTTf) and ifHfl) suggests that F q 

1 /2 

implicitly absorbs an A n factor. When that factor is 
applied to Eqs. O, GSJ) and (gTJ, we expect Eq. H2J 
to follow. This conjecture can be independently checked 
when the centrality dependence of the proton distribu- 
tion at high px becomes available and the existence of 
the scaling §> p (z) can then be examined directly from the 
data, as is done in II. If the conjecture is verified, then 
the data provide empirical evidence for proton being the 
hadronization product of three quarks, rather than other 
mechnisms such as gluon junction pi l22j|. 

One way to check Eq. I|42|) is to calculate the p/w ratio, 
Rp/n, at different centrality bins, where 



R p /tv{PT,N) 



dN,, 



prdp- 



dN K 
PTdp j 



(N) 



(43) 



FIG. 11: Scaling distribution for proton. 



In I this ratio has been calculated for N — iV max . Using 
Eqs. (|41|l and (|42|) we now can calculate it for N < N max . 
There is no data in the pt range where our prediction is 
reliable. However, PHENIX does have preliminary data 
on R p /n for pt < 3.7 GeV/c for two centrality bins, 0- 
5% and 60-91% 15]. They differ by roughly a factor of 3 
with large errors for pt ~ 3 GeV/c. Assuming that the 
ratio of ratios is likely to remain the same at higher pt, 
we can calculate it for comparison, with the definition 

r(Pr) = R p /Apt,N = 350) /R p/n (p T ,N= 15) , (44) 

where N — 15 is taken to correspond to 60-91% central- 
ity. The result is shown in Fig. 12 for y/s = 200 GeV. 
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While the tendency of r(px) to increase at low px is dis- 
turbing, the level of r{px) ~ 3 for p? > 5 GeV/c is in 
rough agreement with the data at pt ~ 3 GeV/c. If 
A P (N) were the same as A n (N), then r(pr) would be 
much lower by a factor of 8.4 and can be ruled out even 
by the preliminary data 4.5]. Thus, until sufficient data 
become available to test directly the scaling formula l|41|) 
for protons, we shall use Eq. i|42|) for the normalization 
factor. 




FIG. 12: The ratio of the p/n ratio at TV = 350 to the same 
at N = 15 for y/s = 200 GeV. 

On the basis of Eq. (|41l) we now can calculate the pt 
distributions, dNp/pxdpx, that are experimentally mea- 
surable. They are shown in Fig. 13 for N = 30, 150, and 
350 for y/s = 200 GeV. It is interesting to note that the 
distributions for N — 150 and 350 differ only slightly in 
the log scale, due undoubtedly to the near cancellation of 
the two opposing properties: the increase of the number 
of hard collisions at higher N and the suppression of high 
Pt protons in larger dense medium. 




FIG. 13: The pr distribution for proton at three values of N 
for = 200 GeV. 



V. SCALING DISTRIBUTION OF KAONS 

For the production of kaons in the recombination 
model we need two inputs: the strange quark distribution 
and the recombination function for a K meson. The for- 
mer has been considered in our study of the strangeness 
enhancement problem in heavy- ion collisions |23j |. The 
latter is given in Ref. [hH ]. 

It is known that the number of A and £ hyperons pro- 
duced in heavy-ion collisions is enhanced by more than a 
factor of 2. Since the enhancement is due mainly to gluon 
conversion and Pauli blocking, the relevant question in 
the context of the parton model is what the percentage 
of the gluon conversion into the strange quarks is. The 
amazing answer found in Ref. j2^| is that it is only 8%. 
Due to the large number of gluons produced in heavy-ion 
collisions that is enough to raise the strange quark den- 
sity by a factor of 2.1 to 2.3 from the intrinsic level in a 
free proton, depending on collision energy. 

To be precise, let us use the notation where the sym- 
bols £ and s denote the total number of light and strange 
quarks, respectively, £ being u 4- d, in contrast to q being 
u or d. Let 4,4, s s and g denote the numbers of valence 
quarks, light sea quarks, strange quarks and gluons, re- 
spectively, before the gluons are converted to ££ and ss 
for recombination to form hadrons. Since only the ra- 
tios of these numbers will be relevant below, they will be 
given modulo a common multiplicative factor that need 
not be specified. The numbers for RHIC at 130 GeV are 
M 



4 = 0.30 , 4 = 0.37 
s s = 0.18 , g = 2.46 . 



(45) 



It is also found that the fraction of gluon conversion to 
the strange sector is 7 = 0.08, i.e., 



19 



4 = (1 - 7)5 , 



(46) 



where the subscript c denotes converted quarks. The net 
strange to light quark ratio after conversion is then 



0.128 



£ 4 4~ 4 4~ 4 

Setting u = d for simplicity, we have £ — 2q and 
s/q = 0.256 . 



(47) 



(48) 



Since s c — jg — 0.197, the strange quark multiplicity 
originally at s s is more than doubled by gluon conversion. 

The above consideration is at the quark level. How 
that translates to hadron abundance must take into ac- 
count hyperon production in addition to kaon produc- 
tion, since the effects of associated production cannot be 
ignored. The problem of partitioning the total strange 
quark numbers into various channels of strange hadrons 
competing for those quarks has been treated in Ref. [2^| . 
It is found there that the fraction k of s quark forming 
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K and k of s antiquark forming K can be deduced from 
the data. That is, defining 



K = KS , K = KS 



(49) 



where K = K~ + K and K = K+ + K° denote the 
numbers of K mesons of various types, one has at 130 
GeV collision energy 



k = 0.628 , R = 0.713 



(50) 



For the purpose of calculating K + /tt + ratio, let us use 
s' to denote the number of s quarks to recombine with u 
quark to form K + , and we have 



- s > = = 0.091^ • 



(51) 



We shall assume that the z distribution of the s' quark 
is the same, apart from normalization, as that of the u 
quark, so we get 



F s >(z)=0.091F q {z) 



(52) 



In the recombination model we expect the produced 
K + to have also a scaling distribution 

®k(z) = / dzidz 2 z 1 z 2 F q (z 1 )F sl (z 2 )Rk(z 1 , z 2 , z) . 



. (53) 

The recombination function for the K meson is similar 
to that of the pion given in I 

R k (z 1 ,z 2 ,z)=R° k z- 2 Gk (-,-), (54) 

V z z / 

where Gk(£i, £2) is the valon distribution in the K meson 


Gjc(£i,60= Sa:^|««i+6-1) (55) 

with gx = B(a + 1,6 + 1) . The parameters a and b 
are determined from the analysis of Kp collisions and are 
found to be a = 1 and b — 2 (see the second paper in 
Ref. [HI)- Since K + is the only state in the pseudoscalar 
octet that has us content, we have B? K = 1 and thus 



Rk(z 1 ,z 2 ,z) 



1 Z\Zo . (Z\ , z 2 



5(2,3) z 5 \ z z 



5 — + — - 1 . (56) 



We now can calculate 3>k(z) using Eqs. {TSJ, J52J|, (|S^|l . 
and |JSSJ. Since for $k{ z ) w e have, as in Eq. (|41|l . 



=A Jf (JV)^( a ,JV)^-( S , AO , (57) 

where Ar: = and = , we obtain for the K + /ir + 
ratio 



Rk/k(Pt,N) 



dN K I dN„ $k(z(pt,N)) 



p T dpT/ prdpr ® n (z(p T ,N)) 

(58) 



where ir + is taken to be the same as ir°. 

In Fig. 14 we show R K/7r (p T ) for N = 350 and 200 and 
y/s = 130 GeV. Since the determination of k and R in Eq. 
(15011 is by use of the data on particle ratios, which are not 
reliable for non-central collisions, we have no confidence 
in the strangeness enhancement factor deduced when N 
is low. For N > 200 Fig. 14 shows that the K/ir ratio 
is not sensitive to centrality, and only mildly dependent 
on pt- This result awaits direct check by experimental 
data. 
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FIG. 14: The K/tv ratio at two values of N for ^/i = 130 
GeV. 



VI. CONCLUSION 

From the scaling distribution <fr(z) obtained from the 
data on ir° production, we have derived many quantities 
by use of the recombination model. They are the distri- 
butions of quarks, protons, and kaons, and their respec- 
tive dependences on centrality. All that is made possible 
by the discovery of the universal function <&{z) valid for 
all centrality and collision energy. The scaling variable z 
that unites the pt dependences for all N and y/s quanti- 
fies the difficulty of producing transverse motion, and for 
that interpretation it can be termed transversity. A par- 
ticle produced at a particular p? at high yfs has a lower 
transversity than that of the same particle produced at 
the same pr but at a lower yfs, because it is easier in 
the former case that involves a smaller momentum frac- 
tion. Similarly, it is easier to produce a particle at a given 
high pt when N is low than to do the same at a higher 
N because there is less degradation of the transverse mo- 
mentum at lower N . Thus the former particle has lower 
transversity than the latter. With this way of viewing 
transverse motion we see that it is the transversity of 
the produced particles that has the universal property at 
whatever centrality and collision energy. 

The quantitative results that we have obtained have 
their limitations due to the assumptions that we have 
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made. For example, for the proton z distribution we 
have assumed that Fq(z) / F q (z) is a constant, which has 
some phenomenological support in that the observed p/p 
ratio is roughly constant in pt', however, that ratio fails 
to maintain constancy at very peripheral collisions. Thus 
our result is not likely to be valid when N is very low. 
Similarly, whether F s i(z)/F q (z) is a constant and over 
what range of N, if it is, are not known. The central- 
ity dependence of the strangeness fraction of gluon con- 
version is also unknown. Our result on K/n ratio can 
therefore only be regarded as preliminary, pending ex- 
perimental guidance to improve our simplifying assump- 
tions. Despite these uncertainties for non-central colli- 
sions, the recombination model has enabled us to calcu- 
late the scaling behavior of the quark, proton, and kaon 
distributions, the latter two of which are subject to direct 
experimental test. 

For non-central collisions we have only calculated the 
Pt distribution, averaged over the azimuthal angle <j). 
Clearly, the dependence on (f> is important as it contains 
dynamical information. It will be very interesting to in- 
vestigate whether centrality scaling persists in restricted 
<f> bins, and if it does, how the scaling curves depend on 
(j). On the basis of the universality in transversity distri- 
bution, we expect that such ^-dependent scaling curves 
can be put into an overall (^-independent scaling curve 
upon rcscaling. At the price of lower statistics this can 
be checked by appropriate analysis of the ir° production 
data. 

Ultimately, the important issue to focus on is the impli- 



cation of the existence of the scaling behavior on the pos- 
sible formation of quark-gluon plasma. At this stage of 
our understanding a conservative statement that can be 
made is that the discovery of scaling violation might pro- 
vide a strong hint for a drastic change of dynamics, possi- 
bly associated with a phase transition. Without waiting 
for LHC to enlighten us with that possibility, a more ur- 
gent issue to settle is why, if the quark-gluon plasma has 
been created already at existing collision energies, the 
change of the nature of the dense medium does not affect 
the scaling behavior that we have found to be universal 
between 17 and 200 GeV. Either the scaling behavior is 
insensitive to the change, or the change has already oc- 
curred at collision energy less than 17 GeV. If both of 
these alternatives are incorrect, then the only way out 
is that quark-gluon plasma has not yet been created at 
RHIC. These are the unintended, but remarkable, con- 
sequences of the scaling behavior. It is thus paramount 
to understand whether a phase transition can lead to a 
violation of the scaling behavior found here. 
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